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Abstract 

We extend to the special case of Bhabha scattering the " Z-peak subtracted" 
representation previously applied to e + e~ —►//(/ ^ e). This allows us to 
analyze the process at any energy while imposing in an automatic way the 
constraints set by high precision measurements at the Z peak. The procedure 
turns out to be particularly convenient for the search of a certain class of 
new physics effects at variable energy. A few examples are considered and 
the information thus obtained is combined with the corresponding one that 
can be derived from the other e + e~ annihilation processes, both at present 
(LEP2) and at future colliders. This shows that the role of Bhabha scattering 
in this respect can be quite relevant. 

PACS numbers: 12.15.-y, 12.60.-i, 13.10+q 
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I. INTRODUCTION. 



A convenient way of searching for virtual new physics effects in e + e~ annihilation pro- 
cesses has been recently proposed; it consists of the use of the so-called " Z-peak subtracted 
representation" that allows to take automatically into account the severe constraints im- 
posed by the high precision measurements performed at the Z peak by LEP1 and SLC [|TJ. 
This is achieved by choosing as "theoretical inputs" the measured values of the partial Z 
widths Tf and of the effective weak angle sin 2 9 e s together with a(0). For each e + e" — > ff 
(/ 7^ e) annihilation process, all the one-loop standard model (SM) or new physics (NP) 
effects are described by four functions of the energy (^fq 2 ) and of the scattering angle 8 
which are subtracted at q 2 = M\ or at q 2 = in order to take the inputs into account. We 
defer to ref. || for further details. 

This procedure is especially suitable for the search of NP effects which grow with the 
energy. Various applications were made for Supersymmetry , Technicolour , anomalous 
gauge couplings, higher Z' bosons, 4-fermion contact terms and extra-dimensions The 
description was shown |J to be particularly useful when NP is characterized by an effective 
scale A which is much higher than the actual energy range of y/q 2 ~. In the particular case of 
universal 9 independent effects all the information on NP can be then conveyed into three 
constants called, 5z, s , 7 , that can be viewed as the generalization, beyond the Z peak, of the 
T, S 0or 61,3 10 description. 

This representation was not yet applied to Bhabha scattering because of the complication 
generated by the presence of t-channel photon and Z exchanges. The purpose of the present 
paper is to fill this lack. We shall show here in Section 2 that the whole Z-peak subtracted 
formalism can be extended to Bhabha scattering in a very natural way. This will allow to use 
this process, together with the other e + e~ — ► // ones, in order to improve the constraints 
on possible NP contributions, as illustrated in Section 3. In Section 4, we shall show in the 
numerical applications that the gain thus obtained can be substancial. 



II. THE Z-PEAK SUBTRACTED REPRESENTATION 



We shall first summarize the results of the Z-peak subtracted representation, described 
in previous dedicated papers 0, by writing the general e + e~ — > ff (/ ^ e) scattering 
amplitude at one loop as the sum of an effective photon and an effective Z amplitude with 
couplings gyj(q 2 , 9), QyAtf 0), g\j{q 2 -, 9), where the index j denotes either the initial electron 
(j = e) or the final fermion (j = / ^ e). 



<T 1 q 2 - M| + iM z Tz 

^(e + M#eV,0)-si?(^^^ (2-1) 

The aforementioned inputs are taken into account by imposing that the total amplitude 
takes the required value at q 2 = and at q 2 = M|. Following the method of 0, this 
amounts to use a subtraction procedure which allows to write: 

It 



g Ve (q 2 , 9) = V^O) Q e [l + -A a , ef (q 
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gl f {q\ 9) = ^na(0) Q f [l + ^A a , ef (q 2 , 9)} 
9le{q\0)=g\ f { q \9)=U 

9 Z y e = 4 he B.[l - \ReM\ 9) - ^\Q f \V e f(q 2 , 9)} 
9$ f (q 2 ,0) = 7? h f v f [l-\ReM\9)-^-\Q } \V e z p{q\9)} 



v f 



g Z Af {q\9)=4 f hf[l-\Ref{q\ 
with the Z-peak inputs 



(2.2) 



and 



i r 487rr,- n i 

^ Wi-if (2 - 3) 



Vj = 1 - 4|Q^ < 2 ' 4 ) 



where s 2 = 1 — c 2 is the weak effective angle measured through the forward-backward or 
polarization asymmetries in the final channel j, s e = s M = s r and iV, is the colour factor 
with QCD corrections at Z-peak. 

The quantities A^ ef (q 2 , 9), R ef (q 2 , 9), V?f(q 2 , 9), V e z f(q 2 , 9) contain all the g 2 , 9 depen- 
dent parts of the scattering amplitude due to SM or NP at one-loop. In our approach, 
they consist of certain finite combinations of self-energies, vertices and boxes that are 
automatically gauge independent. 

For an additional four fermion amplitude 

v(e + )^Hq 2 , 9) - b(q 2 , ^) 7 5 ] M (e-).«(/) 7 ,[c(g 2 , 9) - d(q 2 , 9)^}v(f) (2.5) 

where a, b, c and d are 0(a) quantities, one easily obtains the corresponding projections on 
the photon and Z Lorentz structures: 

2 m _ Ja(q 2 ,9)-b(q 2 1 9)d e }[c(q 2 ,0)-d(q 2 1 9)v f } 



A a , ef (q 2 ,9) = q 



e 2 Q e Qf 



V e } Z (q 2 ,9) = -(q 2 -M 2 



e 2 hehf 

a(q 2 1 9)-b(q 2 ,9)d e }2S e d e d(q 2 1 9) 



e 2 Q e I. 



3/ 



9) = ~(q 2 - M 2 J C ^ 9) - d ^^ (2.6) 

We have given in ref. and in the Appendix of ref. M the expression of the general 
polarized e + e~ — ► // differential cross section in terms of these four functions. From this 
one obtains, for example, the integrated cross sections a y and the asymmetries A FB j, A LR j. 
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To generalize our approach to the case / = e, Bhabha scattering, it is convenient to 
write the scattering amplitude at one loop as the sum of two (s-channel and t-channel) 
components: 

A ee = Mq 2 ,9)+A t (q 2 ,9) (2.7) 

The procedure that we have illustrated in the e + e~ — > ff (f ^ e) case applies directly 
to the s-channel part of the Bhabha amplitude. In this case we can drop the index ef, as 
e = f and we have V lZ (q 2 ,9) = V 21 ^ 2 ^) = V(q 2 ,9), so that we only deal with three 
independent functions A a (q 2 ,9), R(q 2 ,9) and V(q 2 ,9). 

It is now straightforward to check that the same procedure can be applied step by step 
to the t-channel component. In full generality the latter can be written as 

br,( c +\rv»ni7)(„2 »W^2 



A t (q 2 ,9) =^(e + )7^(9 2 ^)^(e + ).«(e-)7^17/(r^)«(e-) + 



with the t-channel effective couplings: 

1— 



(2.8) 



f Ve (q 2 ,9) = ^Ana(0) Q e [l + -A a (q 2 ,9)] 
9vM\ 0) = 7 I he v e [l - -R(q 2 , 9) - -^\Q f \ V (q 2 , 9)] 
g z Ae (q 2 ,9)=4he[l-lR(q\9)} 



(2.9) 



where the new functions A a (q 2 ,9), R(q 2 ,9) and V(q 2 , 9) are obtained from the previous 
s-channel ones (without bar) by (q 2 < ► t) crossing relations: 

A a (q 2 ,9) = [A a (q 2 ,9)](q 2 — t = -^(1 - cos9); cos9 — , 1 + ^) (2.10) 

and analogously for R, V. 

The general expression of the polarized Bhabha differential cross section obtained from 
the sum of the s-channel (|2.1| ) and t-channel (|2.8| ) amplitudes is given in the Appendix, in 
the form: 

da , da 1 , ,. da 2 , , . da p 

l-PP')- -+(1 + PP')- -+(p l -p)- (2.11) 



dcos9 dcos9 dcos9 dcos9 

where P and P' are the initial e~, e + polarizations. 

We shall write the three differential cross sections as the sum of a Born term and a one 
loop contribution, da i = {da i ) Born + (Ax') (1) . 

The Born term is given by 
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<icos# g 2 g 4 t 2 tq 2 

{n( ^ e (u 2 -t 2 + v 2 (u 2 + t 2 ))(q 2 -M 2 ) M 2 , w 2 (g 2 — Mf ) 



W z yL (i + 6 2 )g 2 [(g 2 -M 2 ) 2 + Mlr|] g 2 (t-M 2 ) t[(g 2 -M 2 ) 2 + M 2 r|] 



u 



+ 



3r e . 2r (t 2 + u 2 )(l + £ e 2 ) 2 + 4v 2 (u 2 - 1 2 



t(t-M 2 ) J v aM/ 1 (l+£ 2 ) 2 [(g 2 -M 2 ) 2 + Mfr|] 

u\l+vlf + Avl} f 2(g 2 -M 2 ) | 1 )}] ^ 



(1 + - Mf) V-M 2 ) 2 + M 2 ri] (t-M 2 ) 



v cfcos# 7 * L t 2 aM/(l+£ 2 )t(t-Mf) 

+ W (l+ii e 2 ) 2 (i-M]) 21 1 J 

^ P ]B om = 47raV £ e 3r 1 1 (g 2 -Mf) 1 ] 



v cfcos# 7 q 2 l (l + ti e 2 fw V t 7L [(g 2 -Mf) 2 + M 2 r|] (t - Mf) J 

3r e 2 1 2(q 2 - Mf) 1 

+ [ [(g 2 -M 2 ) 2 + M 2 r|] ( + (t-M 2 ) ) + (t-M 2 ) 2jl ( • } 

and the one loop contributions of the three functions can be written in the condensed way: 

= ^{(t 2 + ^(g 2 , g 2 ) + (u 2 - t 2 )G 2 (g 2 , g 2 ) 
+u 2 [Gi(f , f) + G 2 (t, t) + 2G 1 (g 2 , t) + 2G 2 (g 2 , t)]} (2.15) 

rlrr 2 

{ d^9 ){1) = ™ V[G!(t, t) - G 2 (t, t)] (2.16) 

(-^) (1) = ^^[G 3 (g 2 , g 2 ) + G 3 (g 2 , t) + G 3 (t, g 2 ) + G s (t, t)} (2.17) 
acosu q z 

with 

A a (x) + A a (y) 



Gi(x,y) 



xy 

3r e . v 2 . A a (x) - i2(y) - ^V(y) , A Q (y) - - ^V(x) 



«M/(1+^) 1 x(i/-M 2 ) l/(x-M 2 ) 

3r e ^ R(x) + R(y)+f*$(V(x) + V(y)) 
-) I TZ 5T72T77. iTm J/ 



v aM/ 1 (a;-M 2 )(y-M 2 ) 



3F 

G 2 (x,y) = (— f) 



1 



A a (s) - R{y) + A a (y) - 



3r„ 



aM/(l + i 2 ) 1 x(y-M%) y(x - M 2 Z ) 
4v 2 r + R{y) + *fHV(x) + V{y)) 



a>Mz (1 + ^ 2 ) 2L 



(x-M 2 )(y-M 2 ) 



"]} 



(2.19) 
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G 3 {x,y) = jz ■ .,, {(-^-) 



(1 + ^) «W ar(y-Aff) 

3 r e ^ 2r + J2(y) + ^v(x) + 



■]} (2-20) 



y aM z J 1 (x-M 2 )(y-M 2 ) 

where we use a condensed notation (for x, y corresponding to q 2 or t), A a (q 2 ) meaning 

A a (q 2 ,9) and A a (t) meaning A a (q 2 ,9); similarly for R and V. 

^From the three previous quantities da 1,2 > p we can compute for instance the unpolarized 
angular distribution 

da da 1 da 2 

[2. 21) 



dcosO dcosO dcosO 
the Left-Right polarization asymmetry 

^•») = l5=S>/fes> (2 ' 22) 

and the new (LL+RR) / (LR+RL+LL+RR) polarization asymmetry which arises from the 
typical t-channel scattering amplitude 



III. APPLICATIONS TO SEVERAL NP MODELS 

A. Universal NP with a high scale 

The previous representation eqs.( [2.15"D -( [2.20| ) continues to be valid in the presence of new 
physics (NP) that does not add extra Lorentz structures to those of the SM. In this case, 
one simply decomposes the three general one-loop functions as: 

A Q , R,V= (A Q , R, V) SM + (A a) R, V) NP (3.1) 

and compute the (NP) effects on the various observables, once their contribution to 
(A a , R, V) is specified. 

For a model of new physics that does not satisfy special simplicity requests, the cal- 
culation of virtual effects in the Bhabha scattering is affected by a proliferation of terms 
with respect to the annihilation process e + e~ — > ff, (f ^ e), as one sees immediately from 
inspection of eqs.( |2.l5| )-( j2.20| ). In fact, after ^-integration, one will find in general a set of 



different functions of q 2 that correspond to each power of 9 in the integrand. Each set arises 
from the six original functions A Q , A a , R, R, V, V, which means to double the corre- 
sponding number of the case f^e. Although this can be a purely computational problem, 
it obviously complicates the practical treatment for this process. 



6 



The situation shows a drastic change for those models of new physics that satisfy the 
requests of being, at the same time, universal, independent of the s, t channels scattering 
angle (e.g. only contributing self-energies and/or vertices), and endowed with an intrinsic 
scale A "sufficiently" larger than y/q 1 . In fact, in the Z-peak subtracted approach, one has 
by construction 

A Q (0, 6) = A a (0, 6) = R(M 2 Z , 9) = R(M 2 Z , 6) = V(M 2 Z , 9) = F(Mf , 6) = (3.2) 

For Universal New Physics effects of the previous considered type one can then write the 
following parametrization: 

R UNP {z) = ^=^k Sz] (3.3) 

Z 



t 2 z 



V^(z)=^ J ^-[5 s ] (3.4) 



K NF {z) = w [ 5,] (3.5) 



where z = q 2 , t. 

The quantities 5z, sn will be in general unknown functions of z. For A 2 >> z we can 
reasonably assume that the three functions are smooth. This means that they could be well 
approximated by the coefficient of the lowest power in a q 2 expansion that is Si(0) whenever 
Si(0) 7^ (this will be the case in the two considered examples). In this case, the same three 
parameters will describe the NP effects both on the s and on the t channel observables. 
These parameters are also the same that appear, for the chosen models, in all the remaining 
proccesses e + e~ — > //, (/ ^ e). This fact allows to combine the theoretical analysis of the 
two types of processes without increase of parameters, thus improving the accuracy of the 
conclusions that are reached. 

The NP expression of the functions Gi(x,y) acquires in this case the simple form: 



aM z 1 + vf x-Ml y-M 2 z ' 
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Gs [X ' V) ~ M 2 l + v 2 e i[ aM z )[ y-M 2 z x 1 



y aM z y-M% x - M| 



The three constants 5 Z , 5 S , 5 7 depend on the chosen model and can be easily determined 
in each separate case. To show how this procedure works in practice, we shall provide the 
expressions of the 5, in a couple of specific cases that meet our simplicity requests. With 
this aim, we have considered the following models: 

(1) Anomalous Gauge Couplings (AGC) 
We used the framework of Ref. |H| in which the effective Lagrangian is constructed with 
dimension six operators respecting SU(2) x (7(1) and CP invariance. As shown in Ref. [|J, 
only two parameters (fow an d /zxb) survive in the Z-peak subtracted approach. The explicit 
expression of the UNP contribution to 5 Z , 5 S and <5 7 are 

M 2 c 2 s 2 
5 z = 8na(-f)(^f DW + ^f DB ) 

5 S = S-Ka(—§-)(^f DW - ^db) 
A z s e c e 

M 2 

5, = Sna^ifow + f DB ) (3.9) 
They satisfy the linear constraint: 

5 z - 1 Z 2s2e 5 s + S^ = 0. (3.10) 



(2) Technicolour 

The second considered model was one of Technicolour type with two families of strongly 
coupled resonances (V and A) 0. The typical UNP parameters are the two ratios Fa/ Ma 
and FyjMy where Fa,v and Ma,v are the couplings and the masses (that in this case play 
the role of the new physics scale A >> q 2 ) of the lightest axial and vector resonances. 
The contribution to 5 Z , 5 S and 5 7 are 

ira 9^Ml , My 

b e C e IVly 

M 2 

6, = -Arra^F 2 (3.11) 
Again, we have a linear constraint in the (5 Z , 6 S , <L) space: 



S 



(3i2) 

and the conditions 

5z,s > 5 7 < (3.13) 



B. Non universal examples 

Strictly speaking, our procedure has been motivated by the possibility of investigating 
models of a special universal type, for which the number of parameters to be determined 
can be suitably reduced. But there exist interesting models of NP that, although not of 
universal type, can be nevertheless described by a very restricted number of parameters. In 
these special simple cases our Z-peak subtracted procedure can be applied, without invoking 
any smoothness assumption, using the more general expressions of eqs.( |2.l5D -( ^.2(J ). In what 



follows, we have considered two cases that seem to us particularly relevant. These are: 
(3) Contact terms 

With the idea of compositeness (but it applies to any virtual NP effect with a high intrin- 
sic scale, for example higher vector boson exchanges, satisfying chirality conservation) the 
following interaction 

4-7T - — - — 

+VRL(%r^ R ) (pl-rM) + VLR(^ L m)(H^ f R )} (3-14) 

was first introduced in [|TTJ for any four fermion interaction (ii — > //); kif = | for i = /, 
kij = 1 otherwise; = (1 — 7 5 )/2 ^ R = (1 + 7 5 )/2 \I>; r] ab are phase factors defining the 
chirality structure of the interaction. Various applications have been made for pure chiral 
cases (ij) = LL or RR or LR or RL (keeping only one r\%j = ±1), as well as for mixed cases 
like VV (tjll = Vrr = Vrl = Vlr = ±1), AA = Vrr = ~Vrl = ~Vlr = ±1), VA 
{Vll = ~Vrr = Vrl = ~Vlr = ±1), AV (vll = ~Vrr = ~Vrl = Vlr = ±1); see ref. |TJ 
for a general discussion. 

In the Z-peak subtracted representation, the effect of this interaction on the e + e~ — > // 
(/ 7^ e) observables is obtained through the following expressions: 

2 

K,ef{q\ 9) = ( e 2Q^Q A2 ) [Vll(1 - 0(1 - Vf) + 7^(1 + 0(1 + Vf) 

+ Vrl{1 + v e )(l - v f ) + vlr{1 - 0(1 + v f )] 
K ef{q , V) - -{ or t \2 — n^ LL + Vrr ~ Vrl - Vlr\ 

e -'3e-'3/ A 

9 s r Tr(a 2 — M 2 ) 

V e f(q 2 , 9) = -( ee 2n l \ z) ) [vll(1 - v e ) - vrr(1 + v e ) + vrl(1 + v e ) - Vlr(1 - «.)] 
2s r ir(n 2 — M 2 ) 

V e Z f \q\ 9) = -( e \!? . 2 Z> )[vll(1 - v f ) - vrr(1 + v f ) - Vrl(1 - v f ) + V lr(1 + «/)] (3.15) 



9 



For each choice of chirality structure, like LL, RR, LR, RL, VV ,AA,V A,AV , the effects 
on the differential cross section for two fermion production can be described by a single 
parameter A. 

In the case of Bhabha scattering, the constraint tirl = tjlr applies, so that the above 
expression can be put in the form of eq. ( |3T3| - |3~5| ) , with: 

= ~^ e^A 2 ' Wll + Vrr ~ 2 Vlr\ 

As e c e TcMl. r . . , . . 

5 « = ~( — J WlM 1 ~ u e) - f]RR\ l + v e) + 2v e ri LR \ 

5, = {^)[Vll(1 - v e ) 2 + ^(1 + v e ) 2 + 2 Vlr (1 - t; e 2 )] (3.16) 



(4) Extra dimensions 

Recently, an intense activity has been developed on possible low energy effects of gravi- 
ton exchange. The following matrix element for the 4-fermion process e + e _ — > ff [|l^] is 
predicted: 



— [erf eft fa - Pi)-(p4 - Ps) ~ e^efYf(P2 - Pi)„(p4 - Ps)/J (3.17) 



A 

A 4 " 1 

For this model one finds in the case of Bhabha scattering the following contributions to 
5 7i z,s that, as one sees, are now genuine q 2 , 6 functions: 

,\q 2 M 2 (v 2 e - 2cos6) 



5. 



7 A 4 e 2 

,Xq 2 Ml v l6~s 2 e ~c 2 

*. = (^)<^ (3-18) 

Illustrations will be given in the next Section with the normalization A = ±1. Note that the 
q 2 factor is purely kinematical and a consequence of the higher dimension of the interaction 
Lagrangian. Note also the presence of the term proportional to q 2 cos9 in the s-channel 
photon coefficient 5 7 , which gives a contribution proportional to t + 2q 2 in the t-channel, 
according to eq. fl2.10p . This contribution will turn out to give the largest effect through the 
interference with the standard photon exchange amplitude. 



Our theoretical description of new physics effects is at this point concluded. The final 
Section 4 will be devoted to a detailed numerical analysis of the information that can be 
derived on the involved parameters by the present (LEP2) and future (LC), using both 
Bhabha scattering and all the remaining e + e~ —>■ ff processes. 



IV. NUMERICAL RESULTS 
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A. LEP2 (present and future) 



As a first application of our approach, we have used some of the LEP2 results on fermion 
pair production in order to constrain the set of constants Sz, s ,-y which fully describes the 
effects of general Universal New Physics. In particular, in the same spirit of a previous study 
||, we have considered the following "non Bhabha" observables: a^ T (the cross section for 
H and t pair production), Afb,^,t (the related forward-backward asymmetries) and 05 (the 
cross section for production of quark pairs, for the five light flavours accessible at LEP 
energies). In addition to these observables we have included the unpolarized differential 
Bhabha cross section, measured in intervals of the cosine of the polar angle of the scattered 
electron. For the muon and tau cross sections and asymmetries and for the hadronic cross 
section the combinations of preliminary results of the four LEP experiments with the data 



collected at center of mass energies of 183 GeV and 189 GeV have been used ||14|| . Although 
based on preliminary results, the combined measurements allow to take advantage of the 
whole data sample produced at LEP2 and therefore to benefit of the reduced statistical error 
and of the proper treatment of the various sources of experimental systematic uncertainty. A 
measurement of the differential cross section for Bhabha scattering with acollinearity smaller 
that 10° has been recently performed with a data sample of approximately 180 pb~ l at the 
center of mass energy of 189 GeV ||T5f. The differential cross section is measured in nine 



uniform intervals of the polar angle of the scattered electron, cos# e -, in the range (-0.9, 0.9). 
The precision of the measurement, which is limited primarely by the statistical uncertainty, 
reaches the level of 1% in the interval of most forward scattering angles. 

This measurement, together with the combined results on muon, tau and hadronic ob- 
servables, has been compared to the Standard Model prediction corresponding to the exper- 
imental signal definition. The deviations of the measurement with respect to the Standard 
Model expectations have then been used to measure and constrain the parameters 5z, 8 S 
and S 1 with a \ 2 fit- In the fit procedure the uncertainty on the reference Standard Model 
prediction itself must be taken into account. The theoretical uncertainties on the Standard 
Model predictions for fermion pair production at LEP2 energies are mainly related to the 
estimate of the large QED corrections. In the case of t + t~ and qq production, the dif- 

ferences between predictions of several semianalytic or Monte Carlo calculations fll6|-|18|1 for 
cross sections and asymmetries are smaller than 1% and, therefore, they are negligible with 
respect to the experimental error. In the Bhabha scattering process the different programs 
17,19j20| providing the Standard Model predictions compare each other at the level of 2% 



in the experimental acceptance. Therefore, in our study, we have assigned a 2% uncertainty 
to the reference Standard Model prediction for the Bhabha differential cross section. This 
uncertainty reflects into an error larger than the experimental one in the region of forward 
scattering, which, as will be discussed in the following, is the most sensitive to New Physics 
effects in 5 7 . 

The results of the analysis are shown in Tab. ([1]). As one sees, the addition of Bhabha 
scattering improves, although not spectacularly, the bound on <5 7 which is constrained by the 
data in the forward scattering angle region, where the Bhabha cross section is dominated 
by the t-channel photon exchange contribution. It is interesting to take in mind that , at 
present, the sensitivity to NP effects in the forward Bhabha cross section is spoiled by the 
theoretical uncertainty on the Standard Model theoretical prediction, which dominates over 
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the experimental error. Should this error be reduced, the role of this measurements would 
certainly be more relevant. 

To give a more quantitative meaning to the latter claim, we have simulated a forthcoming 
measurement at 200 GeV with an overall 400 pb -1 luminosity for each experiment, and 
repeated the previous analysis adding these future data to those available at 183, 189 GeV. 
For consistency, we have assumed in all the three sets of data a central value coincident 
with the SM prediction. The errors are those available at 183 and 189 GeV. At 200 GeV we 
considered two scenarios, one with purely statistical errors and one with their combination 
with a 2% theoretical error on Bhabha scattering. The difference between the two cases 
affects mainly the very forward scattering cone and therefore 5 1 . The results of this second 
analysis are shown in Tab. (|T|). As one sees, in agreement with the qualitative expectations, 
the role of Bhabha scattering is now definitely more relevant in the determination of the 
bound for the photonic parameter 8 7 . 

The same two analyses have been performed for the two models involving contact terms 
and extra dimensions. The results are presented in Tables III, IV. For what concerns the 
contact terms, one should first note that the values of the bounds obtained from e + e~ — > 
qq (without Bhabha) depend strongly on the chirality structure. This comes from the 
interference of the contact amplitude with the standard 7, Z exchange amplitude, which is 
larger when the chirality structures of both amplitudes are close to each other. In particular 
the V A bound is found very low because the standard V A amplitude is depressed by the 
small Ze + e~ vector coupling. One then sees that the effect of the Bhabha process on these 
bounds is generally modest. An opposite situation appears for the case of extra dimensions, 
where the bulk of the effect is provided by the "forward" data. As already mentioned in 
Section IIIB(4), the largest effect on the differential cross section comes from the interference 
of the standard photon exchange with the extra dimension term both in the t-channel, 
followed by the mixed (s-channel)- (t-channel) one, the (s-channel)-(s-channel) contribution 
being much smaller. Once again, this is more clearly visible in the analysis that uses the 
future data at 200 GeV, in the case one assumes no theoretical error in the the Bhabha 
component. 

All the numerical results exhibited in Tabs. (||n]) can be represented graphically. We have 
shown in Figs. the planar ellipses that are obtained by projecting onto the three planes 
(5 S , <5 7 ), (5z, S-y) and (5z, S s ) the 95 % C. L. allowed three dimensional region resulting from a 
global fit of all data in terms of the three parameters fe iSj7 . For completeness, we also show 
the results for the two representative AGC and TC models. As a very preliminary comment 
concerning the latter cases we can notice that, although at a rather qualitative level, LEP2 
data apparently do not particularly support the considered TC theoretical proposal that 
would require 5 S > 0. 

B. LC analyses 

This analysis has been performed in a spirit that is very similar to that used for the 
future 200 GeV LEP2 analyses. In other words, we have assumed a set of measurements 
at y^f = 500 GeV whose central values agree with the SM predictions, and postulated a 
purely statistical error corresponding to a high luminosity of 500 fb~ l . We have added to 
the previous LEP2 "non Bhabha" observables the longitudinal polarization asymmetry A^r 
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for lepton production. A discussion of the important role of this observable can be found 
in [21]. Of course, in principle other measurements e.g. for final b or t quarks could be used. 
We have also assumed 9 angular Bhabha measurements for all the three different observables 
(a, a 1 , a p ) defined in Section (2). More precisely, in each bin [9 min , 8 max ] we have considered 

ycos e min da 
°[e mm ,e max ] = J^ dcos6 

Jcos9 max a cos (7 

\ rcos min dfjP 

A LR ,[0 min ,e max ] = / JZ^ dcos6 i ( 4 - 2 ) 

VPninft***] •'•Wflma- U COS U 

1 rcose min da 2 , . . 

A \\,[e min ,e max ] = / -, ~ d cos 1 4.3 

a [e mm fi max \ Jcose max dcos6 

Tables (|V], |VTP contain the numerical results for the universal and not universal models. 
As a general feature, one notices that in the universal case the limits on all the three param- 
eters 5z,s^ are substantially (a factor of two) improved by the use of Bhabha observables. 
The interesting feature is that, in each case, different Bhabha observables play the crucial 
role. In fact, 5z is mostly affected by An (in both angular directions), 5 S (as one expects) 
by Alr ( in the forward cone), and 5 1 by the unpolarized a (again, for very small angles). 
In the considered non universal cases, the effect is, again, not spectacular (although not 
negligible) for the contact terms. Quite on the contrary, there would be a large (a factor 
two) effect in the case of extra dimensions, mostly due to the unpolarized cross section at 
small angles. For this specific model of new physics Bhabha scattering seems therefore to 
represent a fundamental experimental measurement. This statement is well in agreement 



with the results of a recent numerical analysis of LEP2 data [22 . 
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V. CONCLUSION 



In conclusion, we have shown that, for a class of theoretical models of new physics that 
is certainly not empty, the generalization of the Z-peak subtracted approach to the case 
of Bhabha scattering can be simply performed, leading in general to improvements of the 
information that might be obtained. As a general statement, Bhabha scattering appears to 
be always relevant; for models of universal type, polarized and unpolarized observables play 
a crucial role in the determination of the bounds for the different parameters; in other non 
universal interesting cases, like in particular that of extra dimensions, unpolarized Bhabha 
observables appear to play a fundamental role. 
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Appendix A: General form of the polarized Bhabha scattering 
cross section. 



da , .. da 1 , .. da 2 , , . da p 

(1-PP')- -+(1 + PP')- 7:+(P-Ph n ( 5 - X ) 



dcosO dcosO dcosO dcosO 

with 

dcosO ~ levnf 1 {9ve} + q 2 [(q 2 - M 2 ) 2 + M 2 T 2 z f + K ^ e ^ e } + ( " * K ^ Ae } J 

+ [(g2 _ Mf) 1 2 + M|r|] [(t 2 + « 2 )[(^ } ) 2 + (<A?) 2 ] + 4(« 2 - 1 2 )^?^) 2 ] + 2 ^ t {MMf 
• 2m2 ^2) 2 [(^) 2 + aA?) 2 ] + ,,^" M !! 2 , 21 (^) 2 [(^ ) ) 2 + (9%)>) 



q 2 {t - M 2 Y Ve > LKi)Ve ' vyAe ' J t[{q 2 - M 2 ) 2 + M 2 T 

VmiS^ + ( ^ ))2] [(#e))2 + ( ^ ))2] + ^ { sm2) 

4 ( ^ )4 + ^^ ( ^ ))2[( ^ ))2 + ( ^ )2] + (^|F [(( ^ ))2 + ( ^ )2)2 + 4 ^ ) ^ ) ) 2 )} 



(5.2) 



c^ 2 1 



{&^) 4 + T7r^(^ e ) ) 2 [(^ e ) ) 2 - 



dcosO \&itq 2 H 2 ^ Ve ' t{t-M 2 ) 

+^ !F ([(#e ) ) 2 - a*?) 2 ] 2 )} (5.3) 
daP - u2 < 1 K^ASS) 2 ^?^) + («)[(*i?) 2 + (^) 2 D 



dcosfl 4vrg 2 [(g 2 — M|) 2 + Mf r|] v g : 

1 (|7 (7)A2 r -(Z) (Z)^ g 2 ~M|) 

AflVeJ \9Ve 9Ae ) + .r/o s^2\2 , » ,r2 T^2 1 ^Ve J IflVe 5Ue ) 



q\t-Mir^> v^e^y t[(g2 _ M | )2 + M|r 2] 

((^e ) ^ e ) )[(#e ) ) 2 + G^] + («)[G&) a + (9 { Sf\) 



(t - M 2 )[(q 2 - M 2 ) 2 + M 2 r| 



iMf^JS) + 7r4^(#e ) #e ) [(#e ) ) 2 + G^]} (5-4) 



t(t-M|) vyi/e/ (t-M 2 ) 2 
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TABLES 



TABLE I. LEP2. 95% C.L. bounds on <5 S , Z , 7 resulting from a global fit of combined LEP2 
data. In the first column, the observables included in the fit are <7 MiT , Afb,^,t, 05 measured at 183 
and 189 GeV. In the second column, we add the Bhabha unpolarized differential cross section ( in 
nine intervals of the cosine of the scattering angle) at 189 GeV. The last two columns show the 
improvement of the bounds when only the forward (backward) Bhabha scattering measurements 
are included. 





without Bhabha 


with all Bhabha 


forward 


backward 


Sz 


-0.001 ± 0.031 


0.0064 ± 0.028 


0.006 ± 0.03 


0.0011 ± 0.029 




-0.004 ± 0.032 


-0.0087 ± 0.031 


-0.0084 ± 0.032 


-0.0057 ± 0.031 


5y 


-0.0022 ± 0.0083 


0.00019 ± 0.0074 


0.00014 ± 0.0075 


-0.0019 ± 0.0081 



TABLE II. LEP2. 95% C.L. bounds on S s ,z,-y resulting from a global fit of the simulated 
forthcoming combined LEP2 data. The observables included in the fit are u^ T , Afb^t, 05 at 183, 
189, 200 GeV and, in the second column, also the 9 angular measurements of the Bhabha unpo- 
larized cross section at 189 and 200 GeV. For all the observables we assume that the measurement 
are coincident with the Standard Model predictions. The errors are the experimental ones at 183 
and 189 GeV. At 200 GeV, the errors are the statistical ones associated to a 400 pb -1 integrated 
luminosity per experiment with or without (second and third columns) a 2% theoretical error on 
Bhabha scattering. 





without Bhabha with all Bhabha and 2% th. err 


with all Bhabha 


Sz 


0.014 0.012 


0.012 




0.015 0.013 


0.013 


5y 


0.0038 0.0034 


0.0028 
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TABLE III. LEP2. Bounds on non universal New Physics scale at 95% C.L. resulting from a 
global fit of the present combined LEP2 data. We consider the models of Contact Interactions (in 
the eight cases LL, RR, LR, RL, VV, AA, AV, VA) and Extra Dimensions discussed in the paper. 
The data sets used for the fit are the same as in Tab. (H). Since the central value for 1/A is not 
zero, we show the two different bounds for |A| that are obtained in the two cases A < and A > 
respectively. 



without Bhabha with all Bhabha forward backward 

10-9.9 11-9.2 11-9.2 10-9.8 

7.7-12 8.7-10 8.7-10 7.8-12 

6.5-9.2 16-7.8 14-7.3 8.2-9.2 

7.2-15 12-9.7 11-9.5 8.3-13 

13-20 17-16 16-16 13-20 

16- 13 14-14 14-14 15-13 

17- 8.7 17-8.7 17-8.7 17-8.7 
4-3.3 4.2-3.2 4.2-3.2 4-3.3 

0.69-0.75 0.82-2.2 0.8-1.9 0.77-0.9 



TABLE IV. LEP2. Bounds on non universal New Physics scale at 95% C.L. resulting from 
a global fit of the simulated forthcoming combined final LEP2 data. Data sets are chosen as in 
Tab. (0). 





without Bhabha 


with all Bhabha and 2% th. err 


with all Bhabha 




15 


15 


16 


Arr 


13 


14 


15 


Alr 


11 


16 


18 


Arl 


13 


17 


18 


Ayy 


22 


24 


27 


Aaa 


21 


21 


22 


Aav 


16 


16 


16 


Avm 


5.2 


5.2 


5.3 


Aed 


0.89 


1.2 


1.4 
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TABLE V. Linear Collider. 95% C.L. bounds on 5s, z,^ resulting from a global fit of LC 
data assuming that data are taken at 500 GeV with an integrated high luminosity 500 fb^ 1 . The 
"non-Bhabha" observables included in the fit are crj, Afb,i, Alr,i and cr 5 , where I stands for lepton. 
The "Bhabha" observables are the unpolarized angular distribution and the two asymmetries Alr 
and Am defined in the paper. Each of these three observables is assumed to be measured in the same 
9 bins as in the LEP2 analysis. The central values are assumed to coincide with the Standard Model 
predictions and the errors are purely statistical. The first column is obtained without including the 
Bhabha observables in the fit. The second column is with all the Bhabha observables. The next 
three columns are obtained including in the fit only one of the three a, Alr and An. Finally, the 
last two columns are obtained by using the three observables, but only in the forward (backward) 
cone. 





without Bhabha 


with all Bhabha 


a 


A LR 


A \\ 


forw. 


back. 


10 4 


5z 


2.4 


1.3 


1.9 


2.1 


1.7 


1.6 


1.7 


10 4 


S s 


1.4 


0.8 


1.3 


0.85 


1.3 


0.82 


1.3 


10 4 


5 7 


1.1 


0.56 


0.62 


1.1 


0.9 


0.65 


0.82 



TABLE VI. Linear Collider. 95% C.L. bounds on non universal New Physics scales for the 
considered Contact Interactions and Extra Dimensions models. Data sets and columns meaning 
are as m Tab. flv|). Units are TeV. 





without Bhabha 


with all Bhabha 


a 


A L r 


h 


forw. 


back. 


All 


85 


94 


89 


89 


87 


93 


85 


Arr 


84 


92 


87 


88 


85 


92 


84 


Alr 


66 


120 


120 


66 


87 


110 


110 


Arl 


81 


130 


120 


81 


94 


110 


110 


Aw 


120 


150 


150 


120 


120 


150 


140 


Aaa 


110 


140 


130 


110 


120 


120 


130 


Aav 


130 


130 


130 


130 


130 


130 


130 


A V a 


71 


90 


71 


90 


71 


90 


71 


Aed 


3.3 


5.7 


5.7 


3.3 


3.6 


5.6 


4.5 
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FIGURES 



FIG. 1. Two dimensional projections of the 95% C.L. allowed region in the S s ^z,-y space from 
a global fit of LEP2 results on two fermion production. The observables included in the fit are 
ct^t, Afb,h,ti 05 measured at 183 and 189 GeV and the Bhabha unpolarized differential cross 
section (in nine intervals of the cosine of the scattering angle) at 189 GeV. The inner ellipses are 
the projection of the intersection of the 3 dimensional ellipse and the AGC (dashed line) or TC 
(dotted line) constraints. The small cross marks the axes origin, corresponding to the Standard 
Model case. 
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FIG. 2. Two dimensional projections of the 95% C.L. allowed region in the S s ,z,-y space from 
a global fit of LEP2 results on two fermion production. The observables included in the fit are 
o>,t, Afb,h,t, 05 measured at 183, 189, 200 GeV and the 9 angular measurements of the Bhabha 
unpolarized cross section at 189 and 200 GeV. We always assumed the experimental measurements 
to be coincident with the Standard Model predictions. About the errors, we took the available 
actual experimental errors for measurements at 183 and 189 GeV (like in the previous figure) and 
a purely statistical error at 200 GeV under the assumption of an integrated luminosity of 400 pb -1 
for each of the 4 LEP2 experiments. The inner ellipses are the projection of the intersection of the 
3 dimensional ellipse and the AGC (dashed line) or TC (dotted line) constraints. 
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